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Abstract 

Approaches  for  active  control  of  surge  in  axial 
and  centrifugal  compressors  are  based,  for  their 
majority,  on  sensing  them  after  they  occur  and  then 
taking  appropriate  action.  On  the  other  hand,  some 
studies  in  the  past  have  showed  that  precursors  to 
surge  exist  in  some  cases.  But  very  few  studies  have 
actually  made  use  of  precursors  for  active  control.  In 
this  paper  we  show  experimentally  that  surge 
precursors  exist  in  centrifugal  compressors.  Further, 
we  formulate  a novel  control  scheme  in  which  a real 
time  observer  is  used  for  on-line  identification  of 
magnitude  and  frequency  of  dominant  precursor 
waves.  The  observer  outputs  feed  into  a fuzzy  logic 
control  scheme,  and  the  identified  frequency  and 
amplitude  of  the  precursors  are  used  to  actuate  a bleed 
valve  and/or  a fuel  valve  for  active  control  of 
compressor  surge.  Experimental  results  demonstrating 
the  viability  of  the  overall  scheme  are  presented. 

Introduction 

The  design  and  operation  of  jet  engines  are  faced 
with  many  challenges  that  limit  the  operating  range  of 
those  engines.  Among  those  challenges  are 
aerodynamic  phenomena  that  occur  in  the 
compression  system  such  as  rotating  stall  and  surge. 
Rotating  stall  manifests  itself  as  a region  of  severely 
reduced  flow  that  rotates  at  a fraction  of  the 
compressor  rotational  speed  and  causes  a drop  in 
performance.  Surge  is  a pumping  oscillation  that  can 
cause  flameout  and  engine  damage. 

Due  to  the  importance  of  these  phenomena, 
considerable  effort  in  compression  system  studies  has 
been  focused  primarily  on  the  design  and 
implementation  of  active  control  schemes  to  ensure 
stable  operation  of  the  compression  system  over  a 
wide  range  of  operating  conditions1'9.  However,  those 
schemes,  in  their  majority,  are  based  on  sensing  after 
those  phenomena  have  occurred  and  then  taking 
appropriate  control  action.  But  the  severity  of  those 
phenomena  and  their  impact  on  the  engine 
performance  make  it  highly  desirable  to  have  a 


controller  that  is  able  to  avoid  the  operation  of  the 
compressor  in  a region  where  the  compressor  is 
susceptible  to  those  phenomena. 

Some  studies  focused  on  the  identification  of 
precursors  for  surge  and  were  able  to  show  that  such 
precursors  exist  for  axial  compressors  and  may  or  may 
not  exist  for  centrifugal  compressors  depending  on  the 
compressor  under  consideration10'17.  However,  very 
little  use  of  precursors  was  made  in  the  design  and 
implementation  of  controllers,  primarily  because  of 
the  difficulty  associated  with  on-line  identification  of 
such  precursors. 

In  this  paper,  a previously  developed  real  time 
observer18  scheme  is  used  for  on-line  identification  of 
surge  precursors.  Then,  using  the  outputs  from  the  real 
time  observer,  a fuzzy  logic  control  scheme  is 
synthesized  for  control  of  compressor  bleed  and/or 
fuel  flow  to  the  combustor  so  as  to  recover  or  prevent 
surge  occurrence  while  maintaining  the  maximum 
attainable  pressure  ratio.  The  new  control  scheme  is 
implemented  and  evaluated  on  the  centrifugal 
compressor  facility  in  the  School  of  Aerospace 
Engineering  at  the  Georgia  Institute  of  Technology. 

This  paper  is  organized  as  follows:  First,  the 
centrifugal  compressor  facility  is  described  followed 
by  a brief  description  of  the  real  time  observer  scheme 
used  in  this  study.  This  is  followed  by  a description  of 
the  fuzzy  logic  control  scheme  developed  in  this 
study.  Next,  experimental  results  from  the  centrifugal 
compressor  facility  are  presented  to  illustrate  the 
validity  of  the  overall  approach,  followed  by 
conclusions  and  current  focus  of  our  research. 

Centrifugal  Compressor  Facility 

The  primary  purpose  of  the  centrifugal 
compressor  rig  in  the  School  of  Aerospace 
Engineering  at  the  Georgia  Institute  of  Technology  is 
to  provide  a laboratory  facility  for  control-oriented 
studies  on  surge  and  rotating  stall  in  centrifugal 
compressors.  Therefore,  it  is  imperative  that  the  rig 
exhibits  these  phenomena  at  reasonable  time-scales 
and  yet  be  compact  enough  to  be  housed  in  one  of  the 
bays  available  in  the  experimental  facility. 
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The  centrifugal  compressor  rig  includes  an  inlet 
duct,  a single-stage  centrifugal  compressor,  a 
compressor  discharge  duct,  a special  plenum,  a burner, 
an  exhaust  duct  and  a throttle.  Air  enters  the  system 
from  a conical  opening  at  the  tip  of  the  compressor 
inlet  duct  for  smooth  transition  of  flow  from  the  room 
inlet  duct  into  the  compressor.  The  inlet  duct  has  a 
constant  diameter  of  4 inches  (101.6  mm).  It  is  angled 
and  extended  from  the  air  inlet  duct  to  the  compressor 
inlet.  At  the  compressor  inlet,  the  duct  diameter  is 
slightly  reduced  to  fit  the  diameter  of  the  compressor 
impeller.  The  impeller  directs  the  flow  radially 
outward  into  the  compressor  discharge  duct  that  leads 
directly  into  the  plenum.  Heat  is  generated  in  the 
plenum  by  means  of  a propane  burner.  The  flow 
leaves  the  plenum  through  the  plenum  exhaust  duct 
where  it  has  to  pass  through  the  throttle  (butterfly 
valve)  before  it  finally  reaches  the  room  exhaust  duct. 
The  parameters  of  the  compressor  are  given  in 
Table  1. 

The  single-spool  V-5  supercharger  compressor  is 
belt-driven  by  a 40  HP,  3-phase  inverter  heavy-duty 
electric  motor.  It  can  be  operated  at  variable  speeds 
from  approximately  15,000  rpm  to  60,000  rpm 
without  the  need  to  stop  and  make  changes  or 
adjustments.  At  the  maximum  speed  of  60,000  rpm, 
the  compressor  is  capable  of  achieving  a pressure  ratio 
of  approximately  2.25: 1 . 


Parameter 

Value 

Compressor  Inlet  Area 

0.003888  m2 

Compressor  Exit  Area 

0.002150  m2 

Compressor  Inlet  Duct  Area 

0.008107  m2 

Plenum  Inlet  Duct  Area 

0.002027  m2 

Plenum  Volume 

1 m3 

Plenum  Exhaust  Duct  Area 

0.008107  m2 

Compressor  Inlet  Duct  Length  IE- 
shaped) 

0.889  m 

Plenum  Inlet  Duct  Length 

0.3048  m 

Plenum  Diameter 

1.016m 

Plenum  Exhaust  Duct  Length 

0.3048  m 

Compressor  Maximum  Speed 

60,000  rpm 

Compressor  Nominal  Speed 

50,000  rpm 

Table  1.  Centrifugal  compressor  paprameters. 


Data  Acquisition 


Figure  1.  Schematic  of  the  centrifugal  compressor 
facility. 


The  plenum  is  a large  metal  chamber  of  1 m3  in 
volume.  It  is  made  of  aluminum  that  is  3 inches 
(76.2  mm)  thick,  and  is  designed  to  withstand 
pressures  of  up  to  400  Psi  (27.6  bar).  A self-entraining 
burner  is  coupled  to  the  plenum  at  the  bottom.  The 
fuel  to  the  burner  is  controlled  by  a servo  valve.  With 
maximum  fuel  rate,  the  burner  provides  about  1 5 KW 
of  heat  addition.  A commercial  safety  circuit  which 
overrides  the  operator  is  used  for  safe  operation  of  the 
burner.  A quartz  window  provides  visual  access  to  the 
flame. 

The  rig  is  equipped  with  a throttle  actuator  and  a 
fuel  flow  actuator,  which  are  controlled  by  a 
computer.  Both  actuators  have  a bandwidth  of  roughly 
around  10  Hz.  The  fuel  flow  actuator  consists  of  a 
special  flow  control  valve  driven  by  another 
servomotor.  A schematic  of  the  facility  together  with 
the  control  computer  and  the  data  acquisition 
computer  is  given  in  Figure  1 . 

Real  Time  Observer  Scheme 


The  observer  used  in  this  study  was  previously 
developed  at  the  Georgia  Institute  of  Technology  for 
identification  of  combustion  instabilities18.  The 
observer  is  capable  of  real  time  identification  of 
amplitude,  frequency  and  phase  of  one  or  several 
dominant  waves  present  in  a signal.  A brief 
description  of  the  observer  and  its  implementation  for 
the  case  of  identification  of  one  dominant  wave  is 
presented  here.  For  a more  detailed  description,  the 
reader  is  referred  to  Ref.  18. 

The  developed  observer  assumes  that  the 
measured  compressor  pressure  can  be  expressed 
in  the  following  Fourier-like  series19: 

N 

Pm  (0  = X K (0  sin(®»  (00  + 

bn(t)coS(o)„m] 

whose  time  dependent  amplitudes,  i.e.,  an(t)  and 
bn  (t),  and  frequencies  con  (t)  are  as  yet  unknown 
and  need  to  be  determined  by  analyzing  the 
measured  pressure  signal  Pm(t).  Qualitatively,  the 
determination  of  these  unknowns  consists  of  an 
iterative  solution  process  that  is  schematically 
described  in  Figure  2n>.  At  the  start  of  the 
calculations,^  the  frequency  of  the  unknown 
oscillation  (0  is  assumed  and  substituted  into  the 
integrals  in  Figure  2 and  their  lower  limit  of 
integration.  The  resulting  integrals  are  then 

solved  for  the  unknown  coefficients  <7  and  b . 

The  calculated  values  of  a ,b  and  & are  then 
substituted  into  another  relationship  (not  shown 
in  Figure  2)  that  determines^  an  improved  value 
of  the  unknown  frequency  (O  . The  latter  is  then 
substituted  again  into  the  integrals  in  Figure  2 to 

obtain  improved  values  of  Cl  and  b . This 
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calculation  procedure  is  repeated  until  the 

calculated  values  of  a , b and  G>  converge  to 

final  values.  These  values  of  Cl , b and  are 
then  used  to  determine  the  unknown  amplitude, 
frequency  and  phase.  This  solution  procedure 
determines  the  characteristics  of  the  mode  with 
the  largest  amplitude  first. 


Open  Loop  Results 

The  observer  described  in  the  previous  section 
was  implemented  on  the  control  computer  of  the 
centrifugal  compressor  facility.  The  measured  signal 
and  observed  frequency  and  amplitude  of  the 
dominant  wave  for  different  operating  conditions  of 
the  compressor  are  shown  in  Figures  3,  4 and  5.  Figure 
3 shows  results  for  the  case  of  low  compressor 
loading,  i.e.,  for  the  case  of  nearly  full  throttle 
opening.  The  identified  frequency  varies  randomly  in 
the  range  of  1 to  2 KHz,  indicating  a signal  dominated 
by  noise.  As  the  loading  on  the  compressor  is 
increased  by  closing  the  throttle  (see  Figure  4),  the 
identified  frequency  drops  to  around  100  Hz 
indicating  the  presence  of  dominant  waves  in  the 
pressure  measurements. 

As  the  loading  on  the  compressor  is  further  increased, 
the  compressor  surges  (see  Figure  5),  which  is  evident 
from  large  oscillations  in  pressure.  It  is  also  important 
to  note  the  presence  of  cursor  waves  during  part  of 
each  surge  cycle  as  evidenced  by  the  observed 
frequency  and  amplitude. 

In  order  to  verify  the  performance  of  the 
observer,  an  off-line  Fast  Fourier  Transform  (FFT) 
analysis  was  performed  on  part  of  the  input  signal 
indicated  by  two  vertical  lines  in  Figure  5.  The  results 
are  presented  in  Figure  6.  The  first  plot  in  Figure  6 
shows  the  inlet  pressure  variation  during  the  time 
frame  used  for  performing  the  FFT  analysis.  The 
second  plot  in  Figure  6 is  the  time  dependent 
identified  frequency  output  from  the  observer  and  the 
third  one  is  a plot  of  frequency  spectrum  obtained 
from  FFT  analysis  over  the  selected  interval.  A 
comparison  of  the  identified  frequency  output  and  the 
FFT  results  indicates  good  agreement.  Note  the  peak 
around  90  HZ  in  FFT  results  and  the  observer 
identifies  the  frequency  of  the  dominant  wave  to  be 
around  90  HZ  as  well. 
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Figure  3.  Observer  performance  for  the  low  back 
pressure  case  (15,000  rpm) 


Compressor  Speed:  15  KRPM  Throttle  Opening:  10.5547  deg 
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Figure  4.  Observer  performance  for  the  high  back- 
pressure case  (15,000  rpm) 
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Compressor  Speed:  15  KRPM  Throttle  Opening:  7.145  deg 
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Figure  5.  Observer  performance  during  uncontrolled 
surge  (15,000  rpm). 


Compressor  Speed:  15  KRPM  Throttle  Opening:  7.145  deg 
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Figure  6.  Comparison  of  observer  frequency  output 
and  FFT  results  (15,000  RPM) 


Compressor  Speed:  30  KRPM  Throttle  Opening:  1.8323  deg 
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Figure  7.  Observer  performance  during  uncontrolled 
surge  (30,000  rpm). 


The  influence  of  compressor  rpm  on  the 
frequency  of  cursor  waves  was  investigated  by 
repeating  the  open  loop  experiments  with  the 
compressor  rpm  set  at  20,000,  25,000  and  30,000. 
Figure  7 shows  observer  performance  during 
uncontrolled  surge  with  the  compressor  rpm  set  at 
30,000.  Figure  8 is  a comparison  of  the  observer 
frequency  output  and  the  results  from  off-line  FFT 
analysis  performed  using  data  over  a small  time 
interval  (as  indicated  by  the  vertical  lines  in  Figure  7). 
There  is  a good  agreement  between  the  identified 
frequency  from  the  observer  and  the  FFT  results.  Note 
the  small  peak  around  200  Flz  in  the  FFT,  which 
compares  well  with  the  average  value  of  observer 
frequency  output  of  around  200  Flz.  Though  not 
shown,  similar  results  were  also  obtained  from 
experiments  conducted  at  20,000  rpm  and  25,000  rpm. 
From  those  results,  it  was  concluded  that  the  precursor 
frequency  varied  almost  in  proportion  to  the 
compressor  speed  and  this  fact  was  latter  used  in 
closed  loop  experiments. 
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Compressor  Speed:  30  KRPM  Throttle  Opening:  7.8823  deg 


Figure  8.  Comparison  of  observer  frequency  output 
and  FFT  results  (30,000  rpm). 

Fuzzy  Logic  Control  Scheme 

The  control  scheme  uses  a pressure  transducer 
located  at  the  inlet  face  of  the  compressor.  The 
observer  is  continuously  identifying  the  frequency  and 
amplitude  of  dominant  waves  in  the  pressure  signal.  A 
fuzzy  logic  scheme  is  synthesized  to  make  use  of  the 
identified  frequency  and  amplitude  outputs  from  the 
observer  and  apply  appropriate  commands  to  the 
actuators  of  a bleed  valve  or  a fuel  valve.  The 
objective  of  the  controller  is  to  achieve  maximum 
attainable  pressure  from  the  compressor  while 
avoiding  surge.  It  is  very  important  to  note  that  this  is 
achieved  with  very  little  knowledge  of  the  compressor 
characteristic  map. 

The  fuzzy  logic  scheme  consists  of  a fuzzifier,  an 
inference  engine  and  a defuzzifier.  The  inputs  to  the 
fuzzifier  are  the  frequency  distance  I Af| , and  the 
amplitude  difference,  AA.  The  frequency  distance, 

I Af)  , is  the  absolute  value  of  the  difference  between 
the  identified  frequency  and  the  expected  frequency  of 
precursors  which  depends  on  compressor  rpm;  i.e., 

I Afl  = I Identified  Frequency  - Expected  Frequencyl ). 
Likewise,  AA  is  the  difference  between  the  identified 
amplitude  and  a threshold  amplitude.  These  inputs  are 
fuzzifted  using  triangular  membership  functions  with 


two  fuzzy  sets  (“Small’'  and  “Big")  for  each  input  (see 
Figures  9 and  10). 


Figure  9.  Membership  functions  for  identified  |Af|. 


Figure  10.  Membership  functions  for  the  identified 
precursor  amplitude  (AA). 
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Table  2.  Rule  base  used  for  the  Fuzzy  logic 
Controller. 
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The  fuzzified  values  of  the  frequency  distance 
and  amplitude  difference  are  then  passed  through  the 
inference  engine  consisting  of  user  specified  rule-base. 
The  output  consists  of  three  fuzzy  sets,  viz.,  “Negative 
(N)”,  “Approximately  Zero  (AZ)”  and  “Positive  (P)”, 
(see  Table  2). 

The  fuzzy  membership  of  the  output  is  calculated 
using  Sugeno’s  product  inferencing  procedure20  as 
follows: 


XN  = min(iU5(A^)' 


/M 


4/ 


))  + 


mm(ns(AA),nB( 


XAZ  = min  (^5  (^)  ’^b  ( 


A/ 

A/ 


)) 

l 

))  (2) 


xp  = mm(nB(AA),ns 


A/ 


)) 


The  crisp  value  of  control  is  obtained  using  the 
weighted  average  of  the  fuzzy  output  values. 


control  = w..  x.7 
N N 


WAZ  XAZ 


+ wp  XP  (3) 


where  wN,  wAz  and  wP  are  pre-selected  output 
weights. 


Closed  Loop  Results 


The  middle  plot  in  Figure  1 1 shows  that  as  the 
throttle  is  closed,  the  plenum  pressure  increases  until 
the  system  gets  into  surge  at  which  point  the  plenum 
pressure  starts  oscillating.  Note  the  large  decrease  in 
the  plenum  pressure  during  surge.  Eventually,  when 
the  controller  is  turned  on,  the  surge  oscillations 
diminish  and  the  plenum  pressure  is  recovered  to  its 
maximum  value. 

The  lower  trace  in  Figure  1 1 shows  the  throttle 
position  commanded  by  the  user  and  the  actual  throttle 
position  versus  time.  When  the  controller  is  turned  on, 
although  the  user  commands  a more  closed  position 
for  the  throttle,  the  controller  modulates  the  throttle 
position  just  enough  to  get  out  of  surge  and  recover 
the  maximum  plenum  pressure.  The  throttle 
modulation  would  be  synonymous  to  bleed  valve 
modulation  in  a real  engine. 

The  plots  in  Figure  12  are  similar  to  those  of 
Figure  1 1 with  the  only  difference  being  that  in  this 
case,  the  controller  is  turned  on  all  the  time  to 
demonstrate  the  surge  avoidance  capabilities  of  the 
controller.  In  this  case,  we  see  that  although  the 
commanded  position  of  the  throttle  goes  from  a fully 
open  position  to  a fully  closed  position,  the  controller 
does  not  allow  the  actual  position  of  the  throttle  to  go 
lower  than  what  is  necessary  to  keep  the  system  from 
surging  while  maintaining  the  maximum  plenum 
pressure. 


Closed  loop  experiments  were  used  to  evaluate 
the  performance  of  the  surge  controller  described  in 
the  previous  section.  Two  types  of  control  actuation 
were  considered  in  this  study;  viz.,  bleed  valve 
actuation  and  fuel  valve  actuation.  In  the  case  of  bleed 
valve  actuation,  the  exit  throttle  itself  was  used  to 
emulate  a bleed  valve.  We  feel  that  the  fuel  control 
scheme  offers  inherent  advantages  such  as  greater 
control  authority  and  overheating  protection  of  the 
turbine. 

For  each  control  actuation  scheme,  two  sets  of 
experiments  were  conducted.  In  the  first  set,  the 
performance  of  the  controller  for  recovery  from  surge 
was  evaluated.  In  the  second  set  of  experiments,  the 
performance  of  the  controller  in  avoiding  surge  was 
evaluated. 


Throttle  Control 


Figure  1 1 shows  the  time  dependence  of  the  inlet 
pressure  (which  is  used  as  input  to  the  observer), 
relative  plenum  pressure  rise  and  throttle  position  in 
surge  recovery  tests.  The  upper  trace  in  Figure  1 1 
indicates  that  at  large  throttle  openings,  i.e.,  for  low 
compressor  loading,  the  gage  pressure  at  the  inlet  is 
negative  indicating  that  air  is  flowing  into  the 
compressor.  Then,  as  the  throttle  is  closed,  this 
pressure  starts  to  increase  which  is  synonymous  to  a 
reduction  in  the  mass  flow  through  the  compressor, 
and  this  trend  continues  until  the  system  gets  into 
surge.  Note  that  during  parts  of  surge  cycles,  the  inlet 
gage  pressure  is  positive  indicating  flow  reversal. 


Figure  1 1 . Surge  recovery  using  throttle  control. 
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Figure  12.  Surge  avoidance  using  throttle  control. 
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Figure  13.  Surge  recovery  using  fuel  control. 


Surge  Control  Using  Fuel  Valve 


The  same  controller  as  before  was  used  with  a 
fuel  valve  actuation  with  the  only  change  being  that 
the  output  of  the  controller  was  multiplied  by  a 
constant  gain  to  account  for  differences  between 
control  sensitivities  associated  with  the  throttle  valve 
actuation  and  fuel  valve  actuation. 

Figure  13  shows  the  process  of  surge  recovery 
using  the  fuel  flow  control.  The  first  three  traces  are 
similar  to  those  in  Figure  1 1 . The  additional  trace  in 
Figure  13  shows  the  position  of  the  fuel  valve,  which 
is  indicative  of  the  fuel  flow  rate. 

Initially,  the  throttle  is  open,  the  fuel  valve  is 
fully  open  and  the  controller  is  off.  Then,  the  throttle 
is  closed  gradually  until  the  system  gets  into  surge. 
Once  the  fuel  flow  control  is  turned  on,  it  starts  to 
reduce  the  fuel  flow  to  recover  from  surge  while 
maintaining  the  maximum  plenum  pressure. 
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Figure  14.  Surge  avoidance  using  fuel  Control. 


Figure  14  shows  similar  plots  for  the  surge 
avoidance  case  using  the  fuel  valve  actuation.  In  this 
case,  the  controller  is  turned  on  all  the  time.  The 
throttle  is  closed  gradually  until  the  system  is  brought 
to  operate  in  the  surge  region  while  the  fuel  valve  is 
set  at  a constant  value.  Initially,  when  the  system  is 
away  from  surge,  there  is  no  reason  for  the  controller 
to  adjust  the  fuel  flow.  As  the  system  gets  closer  to 
surge,  the  fuel  flow  controller  starts  reducing  the  fuel 
in  order  to  prevent  the  appearance  of  surge  while 
maintaining  the  maximum  plenum  pressure.  In  a real 
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engine,  reducing  the  fuel  flow  when  the  compressor 
gets  closer  to  surge  is  very  desirable  since  this  reduces 
the  temperature,  thus  protecting  the  engine 
components  from  damage. 

Conclusions 

An  observer  previously  developed  at  the  Georgia 
Institute  of  Technology  is  used  to  identify  the 
amplitude  and  frequency  of  the  dynamic  cursors  that 
appear  in  compressors  before  and  during  surge.  The 
amplitude  and  frequency  identified  by  the  observer  are 
used  as  inputs  to  a fuzzy  logic  scheme  for  active  surge 
control.  The  surge  control  scheme  has  been 
implemented  on  the  centrifugal  compressor  facility  in 
the  School  of  Aerospace  Engineering  at  the  Georgia 
Institute  of  Technology  using  both  throttle  valve  (that 
emulates  a bleed  valve)  and  fuel  valve  actuation 
schemes.  Experimental  results  demonstrate  the 
effectiveness  of  the  surge  control  scheme  developed  in 
this  study  for  both  surge  recovery  and  surge 
avoidance. 

Currently  work  is  in  progress  in  providing  a 
detailed  understanding  of  the  dynamic  cursors  and  in 
exploring  the  potential  of  the  active  control  schemes 
developed  in  this  study  for  full  scale  compressors. 
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PAPER  -11,  J.  Prasad 

Question  (F.  Grauer,  Germany) 

Was  the  detector  influenced  by  the  blade  passing  frequencies? 

Reply 

No,  the  observer  could  be  selected  to  operate  over  a range  of  frequencies,  thus  eliminat- 
ing other  frequencies. 

Question  (F.  Grauer,  Germany) 

Do  you  expect  your  fuel  controller  to  be  slow,  as  was  found  from  another  study? 

Reply 

The  bandwidth  of  the  fuel  actuator  we  have  used  is  10  Hz.  We  believe  similar  results 
could  be  expected  from  full  scale  experiments. 


